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^ I Abstract: There has been a significant surge of interest in Horava's model for 3+1 

G^ I dimensional quantum gravity, this model being based on anisotropic scaling at a z = 3 

^ ' Lifshitz point. Horava's model, and its variants, show dramatically improved ultra- 

violet behaviour at the cost of exhibiting violation of Lorentz invariance at ultra- 
^ I high momenta. Following up on our earlier note, Phys. Rev. Lett. 102 (2009) 251601 

5^ \ [arXiv:0904.4464 [hep-th]], we discuss in more detail our variant of Horava's model. In 

contrast to Horava's original model, we abandon "detailed balance" and restore parity 
invariance. We retain, however, Horava's "projectability condition" and explore its 
implications. Under these conditions, we explicitly exhibit the most general model, 
and extract the full classical equations of motion in ADM form. We analyze both 
spin-2 and spin-0 graviton propagators around flat Minkowski space. We furthermore 
analyze the classical evolution of FLRW cosmologies in this model, demonstrating 
that the higher- derivative spatial curvature terms can be used to mimic radiation fluid 
and stiff matter. We conclude with some observations concerning future prospects. 
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1. Introduction 



Hofava's "Lifschitz point gravity" |l|, ||, a field theoretic quantum gravity model 
based on "anisotropic scaling" of the space and time dimensions, has recently at- 
tracted a tremendous amount of attention. Compared to traditional approaches, 
this model exhibits a vastly improved ultra-violet behavior — it is certainly power- 
counting renormalizable, and arguably actually finite ^, ^. Early discussion 
regarding this model includes g |, |, 0, i |, 0, 0, |15|, IH, P, |19 

a M, M, m, m, m 



20, 21 



25, 26, 27, 28 



, 36 1 . In particular, while 



Hofava's specific model of "Lifschitz point gravity" as described in [|l], |], |^ has very 
many desirable formal features, in its original incarnation one appears to be forced 
to accept a non-zero cosmological constant of the wrong sign to be compatible with 
observation |T^, |19|. Additionally, one is forced to accept intrinsic parity violation 
in the purely gravitational sector of the model [|^, ^. 

This naturally leads one to ask (and this is perhaps the key question from the 
quantum-field-theorist's point of view), just how important is "detailed balance"? 
Is it an essential feature of the model, or is it just a simplifying assumption? Are 
there more general models that let you tune the Newton constant and cosmological 
constant independently? Can explicit parity violation in the pure gravity sector be 
eliminated? Below, we shall provide more details concerning a variant of Hofava's 



model that is much better behaved in both these regards |jT9[ , and which appears (at 
least superficially) to be phenomenologically viable. 

Hof ava has also introduced an explicit constraint which he calls the "projectabil- 
ity condition" [Q, From a general relativists' point of view, this is the con- 

dition that a certain part of the space-time metric, the "lapse function", can be 
set globally to unity. Although at first glance this seems a significant constraint, 
the most important key solutions of the vacuum and cosmological Einstein equa- 
tions (the Schwarzschild, Reissner-Nordstrom, Kerr, Kerr-Newman, Friedmann- 
Lemaitre-Robertson-Walker spacetimes) can all be put into this form — at least 
for the physically interesting parts of those spacetimes |19|, E^, BSl. The key question 



(now from the point of view of a general relativist) is this: Is this "projectability 
condition" an essential feature of the model, or is it just a simplifying assumption? 
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A further issue is that Horava's toy model contains a spin-0 scalar graviton in 
addition to the standard spin-2 tensor graviton |l|, H, Phenomenologically, this 
is potentially risky, and might, for instance, run into constraints from the gravity- 
wave-dominated evolution of binary pulsar systems. Should that scalar mode be 
tuned to zero? Is there any symmetry that would protect this? (See for instance the 
discussion in or more general comments within an "emergent gravity" frame- 
work [^, 0.) Moreover, the toy model is purely gravitational, and it will need 
to be investigated carefully just how to embed matter (and eventually the standard 
model of particle physics) within it P, |^, |2^. Because the gravitational field is still 
completely geometrical, albeit with a preferred frame, it might be that there is no in- 
trinsic difficulty in maintaining a universal coupling to the gravitational field. Would 
there be non-zero signals in Eotvos-type experiments? Apart from these questions, 
the community has already begun to look at such things as the possible genera- 
tion of chiral gravitational waves ||^, possible impacts on cosmological solutions and 
perturbations |^, ^, [TD|, |T^, |T^, |2T], |2^, possible modifications of black hole 
physics p, 0, |13|, 0, 123, g9H31], |2], |3|, ig, the question of "absolute time" [|n|]. 



"emergent gravity" [§, , renormahzability [|| |^, |^, , holography p6|] , and 

more... 



Following up on the brief sketch presented in [19|, we discuss in more detail 



some of the phenomenological implications of our extension of Horava's model. In 
contrast to Horava's original model, we abandon "detailed balance" and restore par- 
ity invariance. We retain, however, Horava's "projectability condition" . Under these 
conditions we explicitly exhibit all nine spatial-curvature terms that contribute to the 
Lagrangian. As determined by power-counting, five of these operators are marginal 
(renormalizable) and four are relevant (super-renormalizable). The classical limit of 
this extended model is phenomenologically much better behaved than Horava's orig- 
inal model, as the Newton constant and cosmological constant can be independently 
adjusted to conform to observation. Once the Planck scale and cosmological constant 
have been factored out, the model is described by eight independent dimensionless 
couphngs, one of which is associated with the kinetic energy (and leads to a scalar 
graviton mode), while the other seven are related to the breaking of Lorentz in- 
variance. We demonstrate how the Lorentz-breaking scales are related to the Planck 
scale, but need not be identical to the Planck scale, and sketch the first steps towards 
a detailed confrontation between this quantum gravity model and phenomenology. 

In particular we begin by extracting the full classical equations of motion in ADM 
form, and after linearizing and gauge fixing, use them to analyze both spin-2 tensor 
and spin-0 scalar graviton propagators around fiat Minkowski space, demonstrating 
that these two propagators are sensitive to different subsets of the dimensionless cou- 
pling constants. This analysis also yields the dispersion relations and polarization 
states of classical gravitational radiation in the weak-field regime. Additionally, it 
provides useful insight in the relation between the "projectability condition" and the 
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presence or absence of the spin-0 scalar graviton. We furthermore analyze the classi- 
cal evolution of FLRW cosmologies in this extended model, demonstrating that the 
modified Friedmann equations are sensitive to yet a third subset of the dimension- 
less coupling constants, and that the higher-derivative spatial curvature terms can 
be used to mimic both radiation fluid ("dark radiation") and dark stiff matter. Thus 
different parts of the "potential term" govern distinct aspects of the phenomenology. 
We conclude with some observations concerning future prospects. 



2. The framework: Anisotropic scaling 
2.1 Lapse, shift, and spatial metric 

To explain (our variant of) Hofava's approach [|I|, ^, ^, the basic idea [1^ is to 
write the spacetime metric in ADM form 

ds^ = -N^c^dt^ + gijidx' - N'dt){dx^ - N^dt), (2.1) 

and then, (adopting k as a placeholder symbol for some object with the dimensions 
of momentum), postulate that the engineering dimensions of space and time are 

[dx] = M-i; [dt] = (2.2) 

In condensed-matter language, this is typically referred to as "anisotropic scaling". 
In particle-physics language, one is implicitly introducing a scale Z, with the physical 
dimensions [Z] = [dxY/[dt], and using the theorists' prerogative to adopt units such 
that Z ^ 1. Ultimately we shall interpret this scale Z in terms of the Planck 
scale, and several closely related Lorentz-sjnnmetry breaking scales. If one prefers 
to characterize this scale in terms of a momentum (, then we have Z = C^^^^ c, and 
we see that in order for dimensional analysis to be useful one cannot simultaneously 
set both Z ^ 1 and c ^ 1. (Attempting to set both Z ^ 1 and c — 1 forces 
both dx and dt to be dimensionless, which then renders dimensional analysis utterly 
impotent, and destroys the ability to perform the desired "power-counting" analysis.) 
Consequently in these "theoretician's units" (Z — 1) one must have 

m = H = II = {^r\ (2.3) 

and one is free to additionally choose 

[g,,] = [N] = [1]; [ds] = [^]-\ (2.4) 

To minimize further algebraic manipulations, it is convenient to take the volume 
element to be 

dVd+i = d'x dt; [d\/,+i] = [k]-''-'. (2.5) 
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Note the absence of any factor of c. (This is simply a matter of convenience.) The 
resulting model will, by its very construction, violate Lorentz invariance; the payoff 
however is greatly improved ultraviolet behaviour for the Feynman diagrams ^ 
|, |, coupled with a well-behaved low-energy limit [|l], |, |, g 0]. 



In fact we shall argue that a suitable extension [|19| of the specific model pre- 



sented by Horava in is (at least superficially) phenomenologically viable, and has 
a classical limit that is amenable to analysis in an ADM-like manner. Thus this 
is one of very few quantum gravity models that has any realistic hope of direct 
confrontation with experiment and observation. (When beginning the confrontation 
with phenomenology we will find it useful to go back to the more usual "physical 
units" (c ^ 1) in which Z C^'^.) 

2.2 Extrinsic and intrinsic curvatures 

Like the volume element, the extrinsic curvature is also most conveniently defined to 
not include any explicit factor of c: 

Kij = ^ {-9^J + V.N, + V,iV,} . (2.6) 

Then [A^*] = [da;]/[dt] = in agreement with the previous choices. Furthermore 

For the intrinsic curvature of the spatial slices we have 

[9^,] = [1]; r,k] = H [r',m] = M^ (2.8) 

the key point being 

[Ri'^^] = [VR'^^^] = [k]'^; [V^R^''^] = [k]\ (2.9) 

3. Determining a suitable action functional 

3.1 Kinetic term 

Consider the quantity 

^{K) = SK {{K'^K,, - K') + ^K'} . (3.1) 

(The general relativity kinetic term corresponds to the limit ^ —>■ 0.) Take the kinetic 
action to be 

(3.2) 



Sk = j 3r{K) ^ N d'^x dt. 
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Again, note absence of any factors of c. (This is again purely a matter of convenience 
to keep the dimensional analysis simple.) Then 

[SK] = [gK][^r-'. (3.3) 

Since the kinetic action is (by definition) chosen to be dimensionless, we have 

[9k] = 14'-'^. (3.4) 

Note that the coupling constant gx is dimensionless exactly for 

d = z, (3.5) 

which is exactly the condition that was aimed for in [|^. In a simplified model 
based on scalar field self interactions, this is exactly the condition for well-behaved 



ultraviolet behaviour derived in Q], and also the result obtained in [|T9|. Once we 
have set d = z to make qk dimensionless, then provided gx is positive one can 
without loss of generality re-scale the time and/or space coordinates to set qk 1- 
(A negative qk would ultimately lead to a wrong sign for the Newton constant, and 
qk = is a physically diseased theory that has no kinetic energy terms.) Note that 
there is very little freedom in choosing the kinetic term: -^{K) will be a generic 
feature of any Hofava-like model. 

3.2 Potential term 

Now consider a "potential" term 

Sr = -J y{g) = - j r{9) V9 N d'^x dt, (3.6) 

where 'f'{g) is some scalar built out of the spatial metric and its spatial derivatives. 
Again, note the absence of any factors of c. Then 

[Sy] = [ng)] M-'~^ (3.7) 

whence 

[y{g)] [nr^. (3.8) 

But to keep the kinetic coupling qk dimensionless we needed z ^ d. Therefore 

[y{g)] - [k^. (3.9) 

But "^{g) must be built out of scalar invariants calculable in terms of the Riemann 
tensor and its derivatives, which tells us it must be constructible from objects of the 
form 

{(Riemann)'^, [(VRiemann)]2(Riemann)'^"^ etc...} . (3.10) 

In general, in d+1 dimensions this is a long but finite list. All of these theories should 
be well-behaved as quantum field theories |^, ^]. In particular, since everything 
up to this point is valid for d+ 1 dimensions, these simple observations verify the main 
claims made in |l| regarding 3-1-1 dimensions, and shows that certain key aspects of 
that article nicely generalize to d + 1 dimensions. 
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3.3 Specializing to 3+1 dimensions 



In the specific case d = 3 we have 

[y{g)] - M^ (3.11) 

and so obtain the much shorter specific list: 

|(Riemann)^, [V(Riemann)]^, (Riemann)V^(Riemann), V^(Riemann)|. (3.12) 

But in 3 dimensions the Weyl tensor automatically vanishes, so we can always de- 
compose the Riemann tensor into the Ricci tensor, Ricci scalar, plus the metric. 
Thus we need only consider the much simplified list: 

|(Ricci)^ [V(Ricci)]^ (Ricci) (Ricci), (Ricci)}. (3.13) 

We now consider a model containing all possible terms of this type, eliminating 
redundant terms using: 

• Integration by parts and discarding surface terms. 

• Commutator identities. 



• Bianchi identities. 



• Special relations appropriate to 3 dimensions. 
(Weyl vanishes; properties of Cotton tensor.) 

To keep the calculation tractable, (especially when it comes to integration by parts), 
we impose Hofava's "projectability" condition on the lapse function |jl], [1^. This 
effectively is the demand that the lapse N{t) is a function of t only, not a function 
of position. (In particular, by re-parameterizing the time variable, this means that 
without further loss of generality one can set 1.) Besides simplifying the 

action, Hofava argues that enforcing the projectability condition = N{t) might 
have other merits 0]: Since in his model the action is invariant only under foliation- 



preserving diffeomorphisms (see also section |5.1.4| ), it would not be possible to use 
gauge transformations in order to set = 1 unless one had already enforced N = 
N{t) at the outset. Taking also into account the special nature of time in the theory, 
this could cause technical difficulties in quantization |jl|, |^ . 

It should be remarked that in standard general relativity this projectability con- 
dition can always be enforced locally as a gauge choice. Furthermore for the most 
physically interesting solutions of general relativity it seems that this can always be 
done globally. For instance, for the Schwarzschild and Reissner-Nordstrom space- 
times this "projectability" condition holds globally (for the physically interesting 
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region) in Painleve-GuUstrand coordinates [|1^, |3^, while in the Kerr and Kerr- 
Newman spacetimes this condition holds globally (for the physically interesting r > 



region) in Doran coordinates |I9| , p8| . "Projectability" is also automatic for FLRW 
cosmologies. Thus in this sense the "projectability condition" does not seem to be a 
significant restriction on the physics. 

However, here we wish to enforce projectability at the level of the action, and 
before any functional variation, and there is a price to pay for this: Firstly, not all 
gauges will be accessible to us at the level of the field equations, as will be discussed 
shortly. Secondly, as variation and gauge fixing do not necessarily commute the 
model we are considering might not be the most general model with all possible 
terms of dimension six. Nonetheless, one can expect that such a model will not 
exhibit any qualitative phenomenological deviations from the most general model, 
apart from those clearly related to the projectability condition, which can be easily 
distinguished. Therefore, for the time being we shall take the purely pragmatic 
approach of retaining projectability as a simplifying assumption unless and until we 
are forced to abandon it. The use of this "projectability" condition is a matter of 
some concern and delicacy, as will be discussed later on. 

After a brief calculation, we find that under these conditions there are only five 
independent terms of dimension [k]^: 

RR'jR^i, RjRhR'^i] RV^R, ViRjkV'R^''. (3.14) 

These terms are all marginal (renormalizable) by power counting [|], ^ 0. In 
Hofava's article 0] only a particular linear combination of these five terms is consid- 
ered: 

C'j C\. (3.15) 

The restriction to this (Cotton)^ term follows as a consequence of Hofava's, (to our 
minds), physically unnecessary "detailed balance" condition Furthermore, as we 
shall soon see, this "detailed balance" restriction ultimately makes it difficult to set 
up a phenomenologically viable model based on Hofava's specific proposal [^. 

If we now additionally add all possible lower- dimension terms (relevant operators, 
super-renormalizable by power-counting) we obtain four additional operators: 

\Kf ■ 1; [nf ■ R] [nf : R^; R'^.^ (3.16) 

This now results in a potential with nine terms and nine independent coupling 
constants. In contrast, motivated by his "detailed balance" ansatz, Hofava chooses 
a potential containing six terms (one of which is trivial) with only three independent 
couphng constants, of the form 

■%orava(5') = {h Cottou + Eiusteiu + go metric) ^ (3.17) 
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Note in particular that Hofava's approach includes the cross-term 



(Cotton) X (Einstein) = (Cotton) x (Ricci). (3.18) 

This is a [k]^ term which explicitly violates parity |1|, (and is the only parity- 
violating term in his model). Because we have abandoned "detailed balance" we do 
not need such a term, and we find it convenient to suppress it. It is worthwhile to 
emphasise that when trying to extend or modify Hofava's original model, it is the 
potential term that has most fiexibility: the existence of some y{g) is generic to 
Hofava-like models, but the details may vary from model to model. (For this reason 
many of the computations below are carried out for generic Y{g).) 

3.4 Full classical action 

Assembling all the pieces we now have 

S = j [^{K) - r{g)] ^ N d^x dt, (3.19) 

with 

y{9) = goe+9iCR+ 92 c + 93 e R^,R'' 

+(?4 i?' + 9^ R{R^jR') + g^ R'jRhR\ 

+g7 RV^R + gs ViRjk V'R^\ (3.20) 

where we have introduced suitable factors of C, to ensure the couplings ga are all 
dimensionless. Assuming (71 < 0, we can without loss of generality re-scale the time 
and space coordinates to set both g^ —>■ ^ and gi —1. The Einstein-Hilbert (+ 
cosmological constant) piece of the action is now (still in Z —>■ 1 theoreticians' units) 

5'eh = j {{K"^K,,-K^)+C'R~goC^}^Nd'xdt, (3.21) 

and the "extra" Lorentz- violating terms are controlled by a total of eight dimension- 
less coupling constants (^, (72, . . . , gs) 



Sly = j [iK^-g2 C R' - 93 C R^JR'^ 

—9a R^ — 95 RiRijR"^) — 96 RjR'kR^i 

-gr RV^R - gg ViRjk V'R^''^y/g N d'^x dt. (3.22) 

This is a perfectly reasonable classical Lorentz-violating theory of gravity, which we 
furthermore know has nice ultraviolet behaviour ^, |1^. Even classically, this 
model certainly deserves study in its own right. 
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3.5 Zeroth-order phenomenology: Recovering general relativity 

While the Z ^ \ "theoreticians" units considered above have been most useful for 
power counting purposes, when it comes to phenomenological confrontation with 
observation it is much more useful to adapt more standard "physical" units (c — > 1), 
in which [dx] = [dt]. The transformation to physical units is most easily accomplished 
by setting {d£)z=\ ~^ C~^(di)c=i- In these units the Einstein-Hilbert (+ cosmological 
constant) piece of the action becomes 

SEn = e j {{K^'K,, - K^) + R - goC^} N d'x dt. (3.23) 

(See for instance equation (21.86) of MTW for a comparison with standard 
general relativity.) From this normalization of the Einstein-Hilbert term, we see 
that in these physical (c ^ 1) units 

(167rG'Newton)~' = C'; ^=^5 (^.24) 

so that ( is identified as the Planck scale. The cosmological constant is determined 
by the free parameter qq, and observationally go ~ 10^^^^. In particular, the way we 
have set this up we are free to choose the Newton constant and cosmological constant 
independently (and so to be compatible with observation). In contrast, in the original 
model presented in |jl| , a non-zero Newton constant requires a non-zero cosmological 



constant, of the wrong sign to be compatible with cosmological observations |12, IS 
(In a more realistic model including matter one would have to calculate the vacuum 
energy density and appropriately renormalize the cosmological constant, it would 
then be better to say that the nett renormalized value of go ~ 10~^^^.) 

3.6 First-order phenomenology: Lorentz symmetry breaking 

In "physical" (c ^ 1) units, the "extra" Lorentz- violating terms become 



-g,C' R' - 95C' R{Rr,R'') 

-ge C"' R'jRhR'i - 97 C"' ^V'i? 

-^8 C'^ ^iRjk V'R^^^ d'^x dt. (3.25) 

The extra Lorentz violating terms consist of one kinetic term, and seven higher- 
spatial-curvature terms. The Lorentz violating term in the kinetic energy leads to an 
extra scalar mode for the graviton [Q, with fractional O(^) effects at all momenta. 
(See further discussion below.) Phenomenologically, this behaviour is potentially 
dangerous and should be carefully investigated. In contrast the various Lorentz- 
violating terms in the potential become comparable to the spatial curvature term in 
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the Einstein-Hilbert action only for physical momenta of order 



C{2,3} 



c 




c 



(3.26) 




V^|5'{4,5,6,7,8} 



or higher. Thus the higher- curvature terms are automatically suppressed as we go 
to low curvature (low momentum). 

Note that in this analysis we have now divorced the Planck scale ( from the var- 
ious Lorentz-breaking scales C{2,3,4,5,6,7,8}5 and that we can drive the Lorentz breaking 
scale arbitrarily high by suitable adjustment of the dimensionless couplings g{2,3} 
and (y'{4,5,6,7,8}- It is these pleasant properties that make the model phenomenolog- 
ically viable — at least at a superficial level — and that encourage us to consider 
more detailed confrontation with experiment and observation. Since the ultraviolet 
dominant part of the Lorentz breaking is sixth order in momenta, it neatly evades 
all current bounds on Lorentz symmetry breaking Q At most one might 
hope to get some observational constraints on g2 and/or g^, (which lead to deviations 
from Lorentz invariance at fourth order in spatial momenta), but even those bounds 
would be rather weak. 

4. Classical equations of motion 
4.1 Hamiltonian constraint 

Varying with respect to the lapse N{t) one obtains the Hamiltionian constraint 



The difference here compared to standard general relativity lies in: 

1. The ^ term in the kinetic energy. 

2. The more complicated form of the potential ^((7). 

3. Finally, because of the assumed "projectability" condition on the lapse N{t) 
one cannot derive a super-Hamiltonian constraint, and must remain satisfied 
with this spatially integrated Hamiltonian constraint. 

We emphasise that this Hamiltonian constraint will be generic to all Hofava-like 
models as it really only depends on the "projectabihty" condition and is completely 
independent of the precise form of the potential y{g)- Furthermore, if (and only 
if) one relaxes the "projectability" condition and permits the lapse N{t,x) to be 
an arbitrary function of space-|-time, would one then obtain a super-Hamiltonian 
constraint 



H 




'X = 0. 



(4.1) 




(4.2) 
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4.2 Super-momentum constraint 

Varying witli respect to tlie sliift A^* one obtains tlie super-momentum constraint 

V.TT*^' = 0, (4.3) 

where the super-momentum is 

d[N3r{K)\ 



dg 



[K'^ - Kg'^ + iKg'^] . (4.4) 



The difference here compared to standard general relativity is utterly minimal and 
lies solely in the ^ term. (See for instance equation (21.91) of MTW for a com- 



parison with standard general relativity.) We emphasise that this super- momentum 
constraint will be generic to all Hofava-like models as it only depends on the form 
of the kinetic term ^{K) and is completely independent of the precise form of the 
potential y^g). 

4.3 Dynamical equation 

By varying with respect to gij one now obtains 

— dt{^ n'^) = -2N {{K^^ - KK'^ + ^KK'^] 

\/g 

+ ^^{K) g'^ + (V„iV-) Tx^^ + l^^Txr 



N 6Sr 



y/g Sgi. 



(4.5) 



This is very similar to standard general relativity (see for instance equation (21.115) 
of MTW for a comparison): 



There is a straightforward extra contribution coming from the ^ term in the 
kinetic energy. 

The only real subtlety lies in evaluating the 5Sy/5gij terms. 

We emphasise that most of the features of this dynamical equation will be 
generic to all Hofava-like models. The specific model dependence is confined 
to the 5Sf/5gij term. 

In contrast to standard general relativity there are no terms depending on the 
spatial gradients of the lapse function — this is a side effect of the "projectabil- 
ity" condition. 
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As usual {K'^y^ = K"^"^ grnriK"-^ = K^^ K^^ . Evaluating the 5Sy/5gij term for our 
specific variant of Hofava's model is somewhat tedious, but since in our model we 
know that Sy is the most general action one can build out of the metric using 0, 2, 
4, or 6 derivatives we can (without calculation) deduce that the "forcing term" 

F« = ^ f^. (4.6) 

is the most general symmetric conserved tensor one can build out of the metric and 
0, 2, 4, or 6 derivatives. Writing 

8 

F'^^J29sC'(F,y\ (4.7) 

s=0 

where Ug is an appropriate integer to get the dimensions correct, an explicit calcula- 
tion identifies the following nine individual terms contributing to the overall forcing 
term: 



9o- 1 



91-- R 



92-. 



9s'- RmnR^"' 



{Foh = -\ 9ij- (4.8) 



(FO., = G,^. (4.9) 



{F2)ij = 2RRij - ^R^gij - 2[ViVj - 9ij^'']R- (4-10) 



3 1 

{Fz)i3 — -{RmnR^"')9ij + '^ Rij'^^^ij^ R — ViVjR+3RRij—4:{R )ij — R 9ij. 

(4.11) 

94: R^ ^ 

(F4),,- = SR^Rij - ^R^Qij - 3[ViWj - 9ij^W- (4-12) 

g^i R {RmnR^"^) 

(-^5)^ = Rij{RmnR^^) + 2R{I^)ij — -R{RmnR^^)gij (4-13) 

+ [V\RR,,) + V^Vn{RR"''')9ij - V'ViiRRjk) - V'Vj{RRki)] 

+ [9ij^^ + ^i^j] {RmnR"'^- 
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„ . TDm JDn JDp . 

(Fe),, = ?,{R%, - ]^{R^mR'^r^R''k)gij (4.14) 

+ \ [^\R%3 + ^nVm{R^r'^9^J ' VV^iR^n - "^n^jiR^n] • 

gr. {VRf = {V^R)g^\VnR) ^ 

{F,\, = -2[g,,V^ - V,V,]V'R - 2{V'R)R,, + (V,i?) (VjR) - ^{VRfg,,. 

(4.15) 

gs- (Vpi^^OlV^i?'"") ^ 

(Fg),, = V'R,j + (VpV^V^i?^^)^?,, - VpViV^RPj - V^VjVR^i (4.16) 
-(Vi/2"^„) (V,i?"™) - i(Vpi?^„,)(VP/?'"")(7,, - 2(V''i?,^)(V,i?,p). 

As usual we adopt the notation {R^Y^ = R^^ Rk^ , and similarly {R^y^ = Rk Ri" . 
Furthermore (Vi?)^ = (Vj-R) (V*-R). The first two terms above, (-Fo)ii and (-Fi)ij, 
are exactly what one would expect for the 3+1 decomposition of standard Einstein 
gravity with a cosmo logical constant. The remaining seven forcing terms character- 
ize violations of Lorentz invariance and from a QFT perspective are utterly essential 
for regulating the high-energy behaviour of the Feynman diagrams [|^, ^, ^ . From a 
classical perspective, these are just higher- (spatial)-curvature terms which are sup- 
pressed at low curvature by suitable powers of the relevant Lorentz-breaking scales 
C{2,3,4,5,6,7,8}- The relcvauce of these observations is that the classical limit has now 
been cast into an ADM-like form, suitable, for instance, for detailed numerical in- 
vestigations (and other purposes). 



5. Scalar and tensor graviton propagators around Minkowski 
space 

5.1 Linearizing the equations of motion 

To extract the graviton propagator we will linearize the equations of motion derived 
above. In the classical limit, this calculation will also yield the dispersion relations 
and polarization states for classical weak-field gravitational radiation. For simplicity 
we will set the cosmological constant go = Q for this calculation and use fiat spacetime 
as a background solution. We then have 

V, = 5.,; "AT, = 0; °Ar = l. (5.1) 
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Now consider linearized perturbations 

9ij ^ Sij + € hij] Ni^enf, N^l + en{t), (5.2) 

then 

gij ^S'^ -eh'^ + 0{e^); ^ e n' + 0{e^)- N^l + en{t). (5.3) 

Now expand the extrinsic curvature Kij, the conjugate momentum iTij, the intrinsic 
Ricci curvature Rij and forcing term Fij as a series in e of the form 

X = °X + e^X + 0(e') (5.4) 

In fact Kij, Tiij, Rij, and all vanish at zeroth order and the first non-trivial 
contribution arises at order e. 

5.1.1 Hamiltonian constraint 

Since Kij = 0(e), it follows that ^{K) = O(e^). Furthermore, explicit inspection of 
the potential shows that there is only one non-trivial term 

ng)-eg^CR) + 0{e'). (5.5) 

So the linearized Hamiltonian constraint is 

= gi J ^Rd^x ^0. (5.6) 

This is a rather weak constraint on the integrated Ricci scalar. Note however, that 
if we were to abandon projectability then in the current situation we would have the 
much more restrictive constraint that 

^R = 0. (5.7) 
We shall have more to say on this point later on in the analysis. 

5.1.2 Supermomentum constraint 

Starting from the extrinsic curvature 

^ij = {^ij - ^i^j - ^j^i] + C'(e^), (5.8) 
one calculates the conjugate momentum 

TTij = - diUj - djTii - (1 - -2d-n}^ + O(e^). (5.9) 

Then the supermomentum constraint becomes 

e^d' \^hij - diUj - djUi - (1 - i)5ij{h - 29 • n}| + 0{e^) = 0. (5.10) 
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That is, dropping irrelevant pref actors, 

d' ^hij - diUj - djUi - (1 - 0^ij{h - 25 ■ n}} = 0. (5.11) 

It is convenient to rearrange this as 

d' i^'hij - (1 - O^ijh] = d\ - (1 - 2C)dj{d- n). (5.12) 

We will soon use an appropriate gauge- fixing to eliminate the RHS, but choose to 
defer that step for now. Here and below we shall use as a convenient shorthand 
for the flat-space Laplacian did^. 

5.1.3 Dynamical equation 

We now linearize the dynamical equation. On the LHS 

dt{^ TT^^) = e]^dt [hij - diUj - djm - (1 - OSij{h - 29 • n}} + 0{e^). (5.13) 

In counterpoint, on the RHS the only non-trivial contribution comes from linearizing 
the forcing term ^F^^, all other contributions are O(e^). Consequently 

^dt [kj - diUj - djm - (1 - O^ijih - 29 • n} } = 'Fij, (5.14) 

which it is convenient to rewrite as 

^ \^hij - (1 - OSijh] = \ [d,h^ + - 2(1 - ^)5i^{d • ^)} + ^Fi^. (5.15) 

Before presenting the linearization of the forcing term to explicitly obtain ^F^^ it is 
convenient to first discuss gauge fixing. 

5.1.4 Gauge fixing 

The project ability condition that we have enforced at the level of the action, i.e. the 
requirement N = N(t), is only part of the gauge freedom usually exploited at the 
level of the field equations. The theory under consideration is still invariant under 
coordinate transformations that preserve the preferred foliation: 

t^t + ex°(t) + ...; x'^x' + ex\t, x) + . . . . (5.16) 

(Note that x^i^) depends only on time, not space). Then as per the usual (and quite 
standard) analysis, when working around a fiat background 

gab ^ 9ab + e [daXb + dbXa] + O(e^). (5.17) 
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Extracting the coordinate-transformation-induced changes for the lapse, shift, and 
3-metric (around a flat background) 

N{t) N{t)-ex%t) + ...; (5.18) 
Ni ^ N, + eXi{t,x) + (5.19) 
9ij 9ij + e [d^Xjit, x) + djXiit, x)] + . . . (5.20) 

and so for the linearized perturbations 

n^n- x°(t); 

^ ~\~ Xi ij^i x^ , 
Kj hij + {diXjit,x) + djXiit,x)}; 
h ^h + 2dix\t,x). 

The so-called synchronous gauge consists of choosing x°(^) so that n{t) 0, and 
simultaneously choosing Xi{t,x) so that Ui 0. In this synchronous gauge: 

• From the supermomentum constraint: 

d'{hij-{l-OSijh^ =0- (5.21) 

• From the dynamical equation: 

^{/i.,-(1-0M} = ^F,,. (5.22) 

But even after this synchronous gauge has been adopted there is still a residual gauge 
freedom 

n = 0; ni = 0; (5.23) 

hij hij + diXj{x) + djXi{x)] (5.24) 

h^h + 2dix\x), (5.25) 

where Xi{x) is time-independent. This residual gauge freedom is essential to Horava's 
prescription for separating the perturbation hij onto scalar and tensor modes 

5.1.5 Scalar-tensor decomposition 

One can re-write the supermomentum constraint as 

dt [d' {K, - (1 - i)5,,h}] = 0, (5.26) 

implying 

d' {K, - (1 - i)5,,h} = K^{x), (5.27) 
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where Kj{x) is some arbitrary time-independent vector field. Using the residual 
time-independent coordinate transformations Xi{^) we can transform 

d - K ^ d- K + 2^d\d-x), (5.28) 

and suitably choosing Xi{^) these residual gauge transformations can be used to set 
Ki{x) — *• as long as ,^ 7^ (and this is effectively what Hofava does Note that 
even after setting Ki{x) — > we are still left with some limited gauge freedom 

Xi{x) = di^{x), (5.29) 

subject to the condition 

did'^^ix) = 0, (5.30) 

as such transformations leave Ki unchanged. Under this last remaining coordinate 
freedom 

hij hij + 2didj^{x); 

h + 2d'^^{x). (5.31) 

Unfortunately the exceptional case ^ = is of most direct relevance to standard 
general relativity, but let us put that technical issue aside for now. (We will return 
to this point below.) If we adopt Hofava's gauge fixing then the super-momentum 
constraint can effectively be time-integrated to yield 

d' {hi, - (1 - i)5i,h} = 0, (5.32) 

and the only remaining dynamical equation is 

i|/i.,-(l-OM| = %. (5.33) 



Taking the trace 

l-\djh = 'F. (5.34) 

This is the linearized equation of motion for the spin-zero scalar graviton. The 
potential y{g), and consequently the forcing term ^F, is at this stage generic. (An 
explicit computation for our particular model will follow shortly.) To extract the 
spin-two tensor graviton Hofava then defines (the equivalent of) 

Hij = hi, - (1 - O^ijh; = H,, - .^^^^6,,H- H = -{2~ S^h, (5.35) 



where Hij is transverse, and furthermore separates out the transverse traceless piece 

2 v ^^' 92 



H.J = H,, + l f 5,, -^]H, (5.36) 
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implicitly defining Hij. Therefore 



H 

2 



or 



Hi 



(2 - 30 



A _ ^ 



(1 - , ^ 

(2-30 ' ' 



+ (1 - 



So the dynamical equation (|5.33|) takes the form 



H = 2 ^Fi 



(2 - 30 



djdj 

Q2 



h. 



(5.37) 



(5.38) 



(5.39) 



But in view of the dispersion relation for h already derived in equation (|5.34|) , this 
can easily be rewritten as 



H. 



1-3 



so 



2 ( ^F,, - h,, 'F ] + 



did. 



3 \i 



Q2 



vn 



(5.40) 



(5.41) 



This is the linearized equation of motion for the spin-two tensor graviton. The 
potential y{g), and consequently the forcing term ^Fij, is at this stage generic. 
(An explicit computation for our particular model will follow shortly.) It is easy to 
check that the RHS of the above is both transverse and traceless, as is required for 
consistency. The analysis of this section has now decomposed the metric perturbation 
hij onto a spin-0 scalar, h, and a spin-2 tensor Hij, with appropriate equations of 
motion for each. The only remaining task is to calculate the linearized quantities 
-^Fij and ^F for our particular model, using the gauge-fixing to simplify terms as 
much as possible. 

5.1.6 The exceptional = case 

In the exceptional case ^ = the supermomentum constraint after going to syn- 
chronous gauge (and before any residual gauge fixing) is 



implying 



dt[d'{h,^-6,,h}] =0, 



{hij - 6ijh} = Kj{x) 



(5.42) 



(5.43) 



where again Kj{x) is some arbitrary time- independent vector field. Using the residual 
time-independent coordinate transformations Xi{^) we can transform 



(5.44) 
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implying 

^ + d^X[r,,]] d-K-^d-K. (5.45) 

So in this exceptional case we cannot eliminate Ki completely, though we can elim- 
inate K^ijY (To do this choose d'^Xi to be any of the "vector potentials" leading to 
the "magnetic field" — This will specify Xii.^) up to a gauge transformation 
di^{x).) Thus, after this next step of gauge fixing, there is some time independent 
scalar $(a;) such that Ki{x) = di^{x) and so the (time-integrated) supermomentum 
constraint becomes 

d'{h,j-6,,{h + ^)} = 0, (5.46) 

and we still have available two residual gauge transformations coming from the quan- 
tity ^(x): 

h,j ^ hij -2didj^{xy, h^h-2d'^^{x). (5.47) 
The dynamical equation is now 

^{/^.,-M} = %. (5.48) 

Taking the trace 

-h = ^F. (5.49) 

To obtain the transverse-traceless mode suitable for this situation we first define 

H + ^ 

Hij = hij - 6ij{h + $); h,j = Hi^ — Sij; H = -2h - 3$, (5.50) 

where is transverse, and again separate out the transverse traceless piece 

H.. = n.. -^ , 

2 V ' 92 



= H., + ^(h-^]H. (5.51) 



implicitly defining Hij. Therefore 
or 

did. 



h,, = H,, - - ^5,, - {2h + 3$) + 5,,{h + $). (5.53) 

So (remembering that $(x) is time independent) the dynamical equation (|5.48| ) now 
takes the same form as for ^ 7^ 0. 

= 2 ^F,. - i (5., -^)h = 2 ^F., + (5., - ^) h. (5.54) 

But in view of the dispersion relation for h already derived in equation ( |5.49| ), this 
can easily be rewritten as 

H., = 2 ^F,, - (5,, -^)'F = H,, = 2 (^F,, - ^5,, 'F^+^['F]. (5.55) 
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So despite extra technical complications for ,^ = 0, the ultimate set of differential 
equations looks very similar — though we shall still see some subtleties arise after 
explicit linearization of the forcing term. 



5.2 Linecirizing the forcing term 

We can now return to the linearization of the forcing term Fij ± 
general form for the forcing term without yet imposing any gauge conditions and 
specializing to a Minkowski background 



^F,;,. Using the 



^{F5)ij 

{Fe)ij 
\F,),, 

\Fs).j 



-2[d,d,-6,,d^]CRy, 

d'CRij) + ^5ijd'CR)-didjCRy, 







2[did,-Sijd']d'CR); 



d\'R,,) 



Similarly for the trace 



4d\'Ry, 

0; 



\F,) 

\F,) 

\F,) 
\F,) 
\F,) 
\Fr) 



0; 
0; 



4d%'Ry, 



5.56) 

5.57) 

5.58) 

5.59) 
5.60) 
5.61) 
5.62) 

5.63) 

5.64) 

5.65) 

5.66) 

5.67) 
5.68) 
5.69) 
5.70) 

5.71) 



But calculating ^R^ and ^i? is simple. Before any gauge fixing 

^R^-{d''h-hmn,mn)- (5-72) 

Applying Hofava's gauge condition 7^ 0, plus synchronous gauge, plus in addition 
the residual gauge fixing), 



^R 



-^{d'^hij - didjh) - ^{didjh) = -- 



d^Hij + ^{5ijd^ + didj)h 
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(5.73) 



This allows us to assert 



^(Fi)jj = — d^Hij + ^ (derivative terms acting on) h; 
''^(-^2)1^ = C (derivative terms acting on) h; 

(-^4)ij 
^(-^5)ij 

^{F^)ij — —-d^Hij + ^ (derivative terms acting on) h; 



— -d^Hij + ^ (derivative terms acting on) h; 





^ (derivative terms acting on) h; 



(5.74) 

(5.75) 

(5.76) 

(5.77) 
(5.78) 
(5.79) 
(5.80) 

(5.81) 



For the trace we have 







\F,) 




\F,) 




\F,) 


= 0; 


\F,) 


= 0; 


\Fe) 


= 0; 


\F,) 




\Fs) 





In contrast, in the exceptional situation ^ = we obtain 



1 1 



d''Hij--{Sijd^ + didj)^x) 



^R^d'^^x). 



(5.82) 

(5.83) 

(5.84) 

(5.85) 
(5.86) 
(5.87) 
(5.88) 

(5.89) 



(5.90) 



Note in particular that in the ^ = case we see that is "frozen in"; because 
of the interplay between supermomentum constraints and gauge fixing all the time- 
dependence in ^i? has been eliminated. 

5.3 Spin-0 scaleir graviton 

It is most efficient to separate the discussion of the scalar graviton into three separate 
cases. 



- 22 - 



5.3.1 General kinetic term 7^ 0): 

Collecting terms from the above, for the spin-0 scalar graviton, when ^ 7^ the 
linearized equation of motion is simply 

y)h = -^ [laid' + (-4^, - ^^3) C'd' + (4^7 - l9s^ C'd'^ h. 

(5.91) 

This is a sixth-order trans-Bogoliubov dispersion relation |^ for the scalar mode h. 
Note that only some of the potential couplings contribute here, {gi, g2, gs, gr, gs), 
and it is a different set of couplings from what we shall soon see is relevant for the 
tensor graviton {gi, g^, g^), or is relevant for FLRW cosmology {go — ge). Combining 
the Hamiltonian constraint in equation ( ^.6| ) with equation (|5.73|) we can at least 
deduce that 

d''h{t,x) Sx = 0. (5.92) 



In combination with the equation of motion for the spin-0 excitation, see equation 
( ^.91|) , this implies 

[ hd^x = 0. (5.93) 



But then J h d^x is linear in time: 

/MM)d3.^/M0,.)d3.-H</M0,.)d'.. (5.94) 

If we now use the very last bit of coordinate freedom, (the function ^{x)), we can 
only set J h{0,x) d^x —>■ 0, but we cannot choose to fix J h{0,x) d^x. So the total 
volume cannot be fixed: 

jh{t,x)d^x^t Jh{0,x)d^x. (5.95) 

Thus h is given by the solution of equation (|5.91| ) subject to the constraint equa- 
tion ( p:95|) . 



5.3.2 Specific general relativistic kinetic term = 0): 

For ,^ = we need a separate special-case analysis. In view of our computations 
above, one merely has to make the formal replacements ^/i — > — $(x) followed by 
^ — s> to obtain 

h = llgid' + i-Ag, - Ig^ C'd' + ( 4gr - Igs] C'd'] $(x). (5.96) 



2" V 2"7 ' V 2 ^ 

Note that this is not a "wave equation". The quantity $(a;) is by construction 
time-independent, so this says that at each spatial point x the scalar mode h{t, x) is 
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undergoing "constant acceleration" . Indeed equation ( ^.961 ) can readily be solved to 
give 

h{t, x) = h{0, x) + h{0, x) t + i/i(0, x)t^. (5.97) 



According to equation ( |5.47|) , we still have just enough remaining gauge freedom to 
set h{0,x) and get 

h{t,x) = h{0,x)t+^h{0,x)t'^ (5.98) 

(the function ^{x) can now be chosen arbitrarily). Note that h{0,x) is just the 
right hand side of equation (|5.96| ). Using the linearized Hamiltonian constraint in 



equation (p.6|) as above now gives 



J d'^^{x)d^x = 0, (5.99) 

and so 

h{t,x)d^x= I h{0,x)d^x = 0, (5.100) 



which is the only constraint on the functional form of h{0, x). Again, the total volume 
of space is not fixed. 

5.3.3 Non-projectable case: 

Let us briefly examine now what would happen if we would decide to relax the pro- 
jectability condition. In this case the Hamiltonian constraint would become a super- 
Hamiltonian constraint = 0, which would imply, for the ^ 7^ case d'^h{t,x) = 0, 
and for the ^ = case d'^^{x) = 0. 

• For ^ 7^ the constraint d'^h = 0, together with a boundedness condition on 
h{t, x) and suitable fall-off conditions at spatial infinity, then implies 

h{t,x)=0, (5.101) 

so in this particular case the linearized scalar mode can be eliminated.^ 

• In contrast, for ^ = the equation of motion collapses to h = for all x, and 



so 



h{t,x) = h{0,x) + h{0,x)t. (5.102) 



So in the ^ = non-projectable case there is not enough residual gauge freedom 
to set h to zero. We appear to be left with an undesirable linear expansion 
mode. 



similar result has subsequently been derived in a cosmological setting in reference BE 
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Overall, the behaviour of the spin- zero mode is certainly disturbing and definitely de- 
serves fuller investigation. (Note in particular that in a different context the authors 



of reference p3| have also encountered difficulties with the spin- zero mode).^ 



5.4 Spin-2 tensor graviton 

Extracting the equation for the spin-2 tensor graviton requires a little more work. 
From the results assembled for ^Fij above, we can assert 

H^, = - [gid' + gsC'd^ + gsC'd^] H,, + ^X,,. (5.103) 

Here Xij is a tensor that is linear in the scalar h, independent of ^, and whose tensor 
structure arises solely from partial derivatives and the Kronecker delta. Furthermore 
by construction Xij is both transverse and traceless. This tells us 

= {a(92")925,,- + b{d'''')didj}h. (5.104) 

But tracelessness implies 3a + 6 = while transversality implies a + b = 0. Therefore 
a = b = and the tensor Xij = 0. So the spin-2 tensor graviton satisfies 

H^j = - [gid' + g.C'd' + gsC'd'] H,,. (5.105) 

This is again a sixth-order trans-Bogoliubov dispersion relation 0], now sensitive to 
a different combination of the coupling constants. This dispersion relation does not 
depend on ^, and careful analysis of the special case ^ = shows this equation for 
the tensor mode continues to hold even for that otherwise exceptional case. 

5.5 Low- momentum and high-momentum limits 

In the low-momentum limit the spin-2 graviton is naively seen to have both phase 
velocity and group velocity Cgpi^g —gi- But we had already noted that in going 
to physical normalization for the Einstein-Hilbert term we need to scale space and 
time so that gx — ^ +1 and gi —1. So, as expected, the physical phase and 
group velocities of the low-momentim spin 2 graviton are simply c^pj^j 1- (This 
is very strongly reminiscent of the way in which the Bogoliubov dispersion relation 
leads in the low-momentum limit to a phononic branch with finite speed of sound. 



See, for example, |^, ^.) In short, the behaviour of the tensor mode is physically 
reasonable. 



^The fact that for ^ 7^ the scalar mode can be ehminated at the level of linear perturbations 
does not necessarily mean it is absent altogether. Subsequent to the original appearance of this 
article, several papers have discussed this issue — as well as further difficulties that arise beyond 
the linearized level Note also that the linear expansion mode found for the C = 

non-projectable case need not be physical. Indeed, it was again subsequently shown in that it 
is essentially a gauge artifact. 
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In contrast, the low- momentum phase velocity (and group velocity) of the spin-0 
graviton is seen to be 

Cspino = -TT^ W^^- (5.106) 



So at low momentum this spin-0 mode is propagating (hyperbolic) for C, € (0,2/3), 
but is non-propagating (elliptic) for ^ < and ^ > 2/3. Elliptic modes are poten- 
tially dangerous in that they correspond to imaginary frequencies and can lead to 
exponential instabihties, (see for instance @)- Note that whether this mode is 
propagating (hyperbolic) or non-propagating (elliptic) can now depend on the mo- 
mentum. For examples of similar phenomena in the "analogue spacetime" framework, 
see p7| , The situation is actually worse than this: In Hofava's parameterization 

e i-A 



2-3^ 1-3A' 



(5.107) 



so a propagating scalar mode corresponds to a negative kinetic energy ghost mode. 
(See especially (4.56) of [Q.) Since the presence of a scalar ghost depends only 
on the kinetic term ^{K), this feature is likely to be generic to all Horava-like 
models. (In fact, if one takes seriously the supposed renormalization group running 
of Hofava's original gravity model |l[] from aA = l/3[^ = 2/3] conformal coupling in 
the ultraviolet to a A = 1 = 0] Einstein-Hilbert coupling in the infrared, then this 
evolution has to take one through the parameter region where the scalar mode kinetic 
energy is negative.) These observations again strongly suggest that it is desirable to 
eliminate the scalar mode if at all possible. 

Turning to high momentum, the dispersion relation in that regime can be either 
sub-linear or super-linear depending on the signs of the appropriate couplings — in 
the "analogue spacetime" framework such dispersion relations are referred to 
as "subluminal" or "superluminal" respectively. (In that framework such modified 
dispersion relations are typically used as ways of probing the physics of Hawking 
radiation or cosmological particle production [^, |5^, but in the present context 



the (momentum)^ and (momentum)^ terms are to be viewed as fundamental physics.) 



6. Application to FLRW cosmology 

In a cosmological setting the forcing terms Fij simplify tremendously, and provide 
one with a simple modification of standard FLRW cosmology. (As we shall soon see, 
the Friedmann equations pick up a few extra terms from the higher-derivative spatial 
curvature terms in '^ig), and the kinetic part of the Friedmann equations is slightly 
modified by the ^ term in ^/'{K).) The 3+1 line element is {k G {—1, 0, +1}) 

ds^ = -c^de + a{tf ( + [dO^ + sin^ 9 d^^) 1 . (6.1) 
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6.1 First Friedmann equation 

The extrinsic curvature of the spatial shces is 



K,^ ^ g^. K = (6.2) 

a{t) a{t) ' ^ ^ 



whence 

^{K) ^ -(6 - 90 ^. (6.3) 

Since the spatial slices are (by the definition of FLRW spacetime) constant-curvature 
hypersurfaces we have 

2k 6k 
R^, = ^ ^.S R=^y (6.4) 

So in any FLRW spacetime 

yy) - na) = c + ^ + "P^^+f'^'''' + ^^'^^-+^f (6.5) 

(Note that the gi and terms drop out due to the translation invariancc in the 
spatial slices, and that because of its definition A:^"+^ = k and A;^" = k^. ) In an 
idealized exact FLRW spacetime, because of the spatial homogeneity, the spatial 
integral can be dropped from the Hamiltonian constraint, which now simplifies to 

^(g) + y{9) = 0, (6.6) 
leading to the first (vacuum) Friedmann equation (currently in Z = 1 units): 



a^_y{_a) 
^ ~ ~6 



C]-,-^- (6.7) 



This version of the first Friedmann equation will hold in any Hofava-like model in 
that it depends only on the symmetries of FLRW spacetime and the general features 
of the model, such as the structure of the kinetric term ^{K), and is independent of 
the specific details of the potential y{a). For our specific variant of Hofava's model: 

3 A ^oC^ , g,kC' , 2{3g2 + g3)ee , 4{9g, + 3g, + ge)k 
^"2V + ^ + + ■ ^^-^^ 

Going over to c = 1 units, and rearranging somewhat 

3\d\k A , 2{3g2 + gs)C-'k' , 4(9^4 + 3^5 + ge)^^^ 
^"2V ^ + ^"3+ + " ^^'^^ 

So even in the absence of explicit matter fields, the g2 and gs terms can mimic the 
effects of radiation pressure ( "dark radiation" ) , of either sign depending on the sign 
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of the relevant coupling constants. (See also 0.) Similarly, the (74, g^, and terms 
can mimic stiff matter of either sign depending on the sign of spatial curvature and 
the sign of the relevant coupling constants. (Recall that for or stiff matter p = p and 
so Pstiff oc 1/a^.) Note that the stiff-matter mimicking term is absent in Hofava's 
original model p. This is due to the fact that the the only [k]^ term in Hofava's 
original model is (Cotton)^, and the Cotton tensor automatically vanishes for the 
constant curvature spatial slices of a FLRW spacetime. 

We now add a cosmological fluid through a purely pragmatic and phenomeno- 
logical approach: Working within the cosmological hydrodynamical approximation 
we approximate the cosmological stress-energy tensor by two quantities, the density 
p and pressure p, and simply add them to the vacuum equations by demanding the 
correct limit as one approaches classical general relativity. (More fundamentally, we 
could attempt to derive the matter contributions from an action principle, but for 
the present purposes that would be overkill. A model for the matter sector based 
on a simple scalar field is outlined in |^, and the resulting Friedmann equations 
[insofar as there is an overlap] are compatible with our purely pragmatic approach.) 
Under these assumptions the first Friedmann equation becomes 

3 \d\k A , 2i3g2 + gs)C'k' , 4(9^4 + 3^5 + ^?6)C-^A: , p 
'-2y ^ + ^=3+ + + -^' ('-'^^ 

or equivalently 

3 \ a\ k _A 2(3^2 + ^3)r'fc^ 4(9^4 + 3^5 + ^6)C~'fc SrtG^p 



2V a2 a2 3 '3 

(6.11) 

We this see a controlled deviation from standard cosmology, with the deviations from 
standard cosmology being governed (effectively) by three parameters: ^, {3g2 + gs), 
and {9g4 + 3g5 + go). For a generic Hofava-like model coupled to matter we would 

^^^^ 3.\ r(a) SttGnP 



It is worth pointing out that the presence of the second and third term in the right 
hand side of the first Friedmann equation ( |6.11| ) could potentially lead to bouncing 



solutions for suitable values of the parameters. However, the behaviour of matter 
contribution close to the bounce would be critical and further investigation is needed 
(see also [|l3l)- In this regard it is perhaps worthwhile recalling that a bounce requires 



overall violation of the strong energy condition (SEC) in its immediate vicinity 1^2 



53|. 



Before going further, a remark is appropriate concerning the passage from the 
Hamiltonian constraint to equation ( p.lll) : Dropping the integral over space depends 
crucially on exact homogeneity (hence the emphasis on the word "exact" there). 
Therefore, this step will not be possible in more general cosmological contexts, where 
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the universe is not exactly FLRW, but can instead be interpreted as a large number 
of almost FLRW sub-regions. In this case one would have to develop a more complete 
treatment.^ 



6.2 Second Friedmann equation 

The second Friedmann equation comes from the djTiamical equation for the canonical 
momentum ir^^ . For the conjugate momentum in a FLRW universe 

^ _(2 _ 30 - (6.13) 
a 

and so equation ([4. 5]) takes the form 



-a-'dt I a' {2- 30 ^ gA = 2(2 - 30 ^ 9^' + \-^{K) g^^ + F'^ . (6.14) 



This leads to 



Taking the trace 



But from the definition of the forcing term Fj^, its trace satisfies 

So the dynamical equation for the conjugate momentum becomes [Z = 1 units, still 
for vacuum) 

In physical (c — 1) units the only change is 

Adding a (phenomenological) cosmological fluid simply adds a pressure contribution 
3 A a If 3 A .1 d\r(a)a^] , ^ , , 



^Subsequent to the original appearance of this article, Mukohyama Q has advocated the view- 
point that, in such a setting, the integrated Hamiltonian constraint can be re-interpreted as the 
ordinary Friedmann equation with the addition of "dark dust" that is comoving with the preferred 
foliation. 
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We again see a controlled deviation from standard cosmology, with the deviations 
from standard cosmology being governed by ^ and by "^(a) which is at this stage of 
the calculation arbitrary. 

For our specific variant of Hofava's model the explicit forcing terms F^j also 
simplify tremendously in FLRW cosmologies. In particular 

{Fo)ij 
iF2h 

iF5)ij 



1 


(6.21) 


k 






(6.23) 


2A;2 
a\t) 


(6.24) 


108k 
a%t) 


(6.25) 


36k 
a^{t) 


(6.26) 


12k 
a%t) 


(6.27) 


0. 


(6.28) 


0. 


(6.29) 



{F6)ij 
(F7h 

(Fsh 

That is, inserting suitable dimensional factors 

P / goC 9ikC' , 2(3^2 + ^3)A:^C^ , 12(9^4 + 3^5 + ^e)^ 1 

= [—2 ^ + + / ^^-^^^ 

It is then easy to explicitly verify that for our particular variant of Hofava's model 

f 1^,, , ady(a)] 1 d\y(a)a^] 

F, = [-^na) - 9« = -g^^^ 9.. (6.31) 

The deviations from standard cosmology in our model are governed (effectively) by 
three parameters: ^, and the compound parameters {3g2 + g^) and (9gi4 + 35^5 + Qq) 
hiding inside the potential '^{0). 

6.3 Third Priedmann equation 

By eliminating the a?/a^ term between the first and second Priedmann equations 
one sees 
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Going to c ^ 1 units, and adding (phenomenological) cosmological pressure and 
density, 

with y{a) still being arbitrary at this stage. In summary, most of FLRW cosmology 
survives with only minimal changes. There is a small change in the "kinetic" part of 
the Friedmann equations, easily dealt with by inserting a factor of (1 — |^), and in 
the "potential" the higher-spatial-curvature terms mimic various forms of matter. In 
our specific model the y{a) terms mimic radiation pressure ("dark radiation") and 
stiff matter. 



7. Discussion 



The specific extension of Hofava's model that we have outlined above so far only 
considers pure gravity (and a phenomenologically introduced notion of cosmological 
fluid). It is however a very definite proposal with a small number of adjustable 
parameters, (many fewer adjustable parameters than the standard model of particle 
physics), making it worthwhile to put in the additional effort to develop precision 
tests that would confront this model with experimental and observational bounds. 

At this stage we have investigated the graviton propagators (and weak-field grav- 
itational waves) around fiat Minkowski space, both for the spin-0 scalar and spin-2 
tensor gravitons, and demonstrated how they depend on the specific terms in the po- 
tential y{g)- The presence and behaviour of the spin-0 mode appear to be worrying, 
and definitely require further investigation.^ We have also investigated FLRW cos- 
mologies, and seen how the Freidmann equations are modified for generic y{g)- We 
have seen that the Friedmann equations and the graviton propagators are sensitive 
to different subsets of the coupling constants and so can in principle probe different 
parts of the physics. 

The most significant theoretical restriction we have retained is Hofava's "pro- 
jectability" constraint. Our analysis showed that there are a number of reasons why 
it might be useful to relax this constraint. Based on the results of this article we can 
make some suggestions as to what might happen in the absence of a projectability 
constraint: 

• One would have to modify the use of the integration by parts argument when 
discarding surface terms. Some of the terms we discarded when constructing 
our list of all possible terms appearing in ^{g) would have to be retained. (So 
^{g) would then contain significantly more than the five renormalizable and 
four super-renormalizable terms we focused on in the present article.) 

^ After the original upload of this this article we became aware of reference ||5^ , which also seems 
to suggest that serious problems are introduced by the scalar mode. 
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• Without "projectability" the Hamiltonian constraint would become a full super- 
Hamiltonian constraint. 

• The dynamical equation would then pick up additional terms depending on the 
gradient and the Hessian of the lapse. (The forcing term Fij would in general 
be significantly more complicated.) 

• For both 7^ and ^ = the flat-space spin-2 tensor graviton propagator 
we have calculated would not be affected; while the spin-0 scalar graviton 
can either eliminated at the linearized level for ^ 7^ under suitable fall-off 
conditions at spatial infinity, or collapses to a linear expansion/contraction 
mode for ^ = (in the specific gauge used in this article). 

• The Friedmann equations we have extracted would not be affected. (The lapse 

is unity in any FLRW spacetime.) 

Thus generalizing to arbitrary lapse might not be as difficult as initially envisaged [||] , 
and this is a worthwhile option to consider. 

Turning to the future, the most obvious tests of our current variant of Horava's 
model would come from the observational limits on Lorentz violations |4^ , 
but by inspection the model should also fall into the PPN framework, and specifically 



be subject to "preferred frame" effects [54] — this will lead to stringent limits on 



the size of the Lorentz breaking parameter C, arising from solar system physics. Up 
to this stage we have not had to make any specific commitment as to how matter 
(in the form of elementary particles as opposed to cosmological fluid) couples to the 
gravitational field: Because the gravitational field is still completely geometrical, 
albeit with a "preferred frame" , it seems to us that there is no intrinsic difficulty in 
maintaining a "universal" coupling to the gravitational field, and so there is no need 
for violations of the equivalence principle in this class of models — we expect the 
universality of free fall to be maintained and to not see any likelihood for non-zero 
signals in Eotvos-type experiments. 

In conclusion we would argue that this is one of very few quantum gravity models 
that has any realistic hope of direct confrontation with experiment and observation, 
and that it is well worth a very careful look. 
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